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Using the self-consistent Born approximation we calculate Reiter's 
wavefunction for a single hole introduced into the undoped and orbitally 
ordered ground state of the t-J model with ti g orbital degrees of freedom. 
While the number of excitations is similar to the spin t-J model for a given 
J/t, a distinct structure of the calculated wavefunction and its momentum 
dependence is identified suggesting the formation of a novel type of mobile 
polarons. 
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Introduction — Recently, using self-consistent Born approximation (SCBA), 
we showed that a single hole introduced into the undoped ground state of an 
orbital t-J model with ti g orbital degeneracy propagates coherently as a quasi- 
particle [1]. This striking result contradicts the naive expectations which suggest 
that a hole should be trapped in this Ising-like ordered ground state. In fact, 
the motion of a single hole is due to the frequently neglected three-site terms 
and we showed [1] that this new mechanism of hole movement is fundamentally 
different from the coherent hole motion via quantum fluctuations in the standard 
spin t-J model [2], or in the case of e g orbitals [3]. Though, a more detailed 
understanding of this novel mechanism is needed. Hence, instead of considering 
the Green's function of the problem [1] , we investigate the corresponding Reiter's 
wavefunction [4] calculated in the SCBA. 

Orbital polarons — The orbital t-J Hamiltonian, relevant for ferromagnetic 
ab planes with two active t2 g orbitals, zx = b and yz = a, at each site, consists of 



(1) 



2 



three terms [1] H = TL t + Hj + H^s, with 

Ht = - 1 E + 5 K+y + /l - c ') ' ( la ) 

i 

h j = ^Efc^-zW). (lb) 

to) v y 

^ 3s = - i J E (&i-ft"i<A+* + «!-y^(A+y + ^- C ') . ( lc ) 

i 

where a tilde above fermion operators denote the restricted Hilbert space without 
double occupancies, the total on-site density hi — h ia + hn, = a\a i + b\b i , and the 
pseudospin operators Tf = \{fnb — h ia ). In Eq. (lc) we neglected the three-site 
terms which require orbital excitation and therefore do not change the physical 
properties of the system in the low doping regime (cf. discussion in Ref . [1] ) . 

Following Ref. [2] we reduce the t-J model to a polaronic problem, which 
is a physically justified procedure for the AF or AO ordered phases with low 
concentration of added holes [5]. Hence, we: (i) divide the square lattice into two 
sublattices A and B, (ii) rotate pseudospins on the A sublattice, and (iii) introduce 
fermion operators h\ and hard-core boson operators a\ such that 6j = h\at and 
bi = h\(l — a\ai). Then, in the linear spin-wave approximation and having only 
one doped hole in the plane, we obtain the following Fourier-transformed polaronic 
Hamiltonian, H = H t + Hj + H 3s , with 

k,q 

Hj =^oJ2 ( a kA a kA + at B a kB j , (2b) 

k 

^ 3s= E [ e 2/( k ) /l L^kA + ^(k)/i] cB /i kB . (2c) 

k 

Here the sums go over all momenta {k} in the Brillouin zone for the whole lattice* , 
the total number of sites is N, and the orbiton energy is luq = J. 

The vertices and the dispersion relations are equal to (with v = x,y): 

M„(k,q) = 2icos(fc„ - q„), e 1/ (k) = \ Jcos{2k v ). (3) 

Reiter's wavefunction — In the spirit of Ref. [4], let us assume that the 
wavefunction for a hole with quasiparticle (QP) momentum k and initially doped 
into the sublattice A and B, respectively, takes the form: 

|* kA ) = <4 0) (kK A |0) + ^£a«(k,q)4_ qB at JO) 

q 

+ A? E a A ) ( k ,qi j q2)4-q 1 -q 2 A Q! qiA Q! q2i3l ) +••• ' ( 4& ) 

qi,q2 



*One can also perform the sum only in the reduced Brillouin Zone although this requires to 
change — » N/2 everywhere before summations. 
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+ ^ E a B 2) (k,qi,q 2 )^L qi - q2 B«q lB "UlO) + --- • (4b) 

qi,q2 

(0) 

The coefficients a L (k) (with the sublattice index L G {A, B}) are the QP 
spectral weights which follow from the normalization of the wavefunction, whereas 
a^\\s., qi, ...,q n ) for n > are to be determined from the Schrodinger equations 
H\^ kL ) = AkL|*kL)- Substituting Eq. (4a)-(4b) into them yields 

a { J n \]<,{qi}2n)^a { 2 n ^ 1 \h,{qi}2n-i)M^k2 n -i,q2n)GA(^2n,^2 ) )^ ( 5a ) 
a^" _1) (k,{qj 2 n-i) =a^" _2 \k{qj2T l -2)M x (k 2 „-2, q2n-i)G B (k 2 „-i, ^ _1) ), 

(5b) 

a B™\H{q t }2rO=aB™ _1 \k,{qj2n-i)M x (k2„-i,q2„)GB(k2™,M k 2 B ) ), (5c) 
a^" _1) (k,{qj2n-i) =as" _2) (k, {q 4 }2r l -2)Mj / (k 2 „-2, q2n-i)GU(k 2n _i, ^^~ V> ), 

(5d) 

where k„ = k- qi - ... - q„, {qj„ = {qi, q 2 , q„}, ^ = A kL - naj . The 
Green's functions Gl are defined by the self-consistent equations 

G A 1 (k, W )=a>- £ ,(k)-l^GB(k-q,c.-c )M2(k,q) , (6a) 
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G B 1 (k^)=a>- £;c (k)-l^G A (k-q,c.-c.o)M2(k,q) , (6b) 

q 

and the QP energy has to be equal 

A kA = e,(k) + l^G B (k-q, M «)A/4 2 (k,q) , (7a) 
q 

A kS = e x (k) + 1 G ^( k - <1. Mk'^M^k, q) . (7b) 
q 

Note, that in order to obtain equations for the coefficients of the Reiter's 
wavefunction a^f 1 Eqs. (5) we adopted the following contraction procedure: we 
neglected all terms which would correspond to the annihilation of orbitons in the 
n-th step if the orbitons were created earlier than in the (n — l)-th step. This 
procedure resembles the non-crossing approximation while calculating Green's 
function using the diagrammatic technique. In fact, one may wonder why we 
still have to adopt such approximation since the closed loops which correspond 
to the crossing diagrams are anyway prohibited in the orbital t—J model under 
consideration [1]. The answer to this puzzle is the following: to conclude that the 
crossing diagrams are unphysical we not only need to look at the structure of the 
Hamiltonian Eq. (2) but also at the processes allowed by this Hamiltonian on the 
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square lattice. Thus, we need some extra information about the lattice to exclude 
the crossing diagrams, and consequently we need to introduce the contraction 
procedure by hand. 

Consequences — Eqs. (5) together with Eqs. (6) and (7) can be easily solved 

(n) 

numerically. While it is impossible to calculate the coefficients a y L for all n, one 
may ask whether there exists such an m that all the coefficients with n > m 
are so small that they could be safely neglected. This question is also of physical 
importance since knowing m would mean that the wavefunction of the doped hole 
can be approximated by a superposition of the wavefunction of a free hole and m 
wavefunctions of the hole dressed with m orbitons. 

The easiest way to answer the above question is to calculate the norm 

N^ A = (* kA |*kA) = |a( 1 0) (k)| 2 {l + l2|M !B (k,q 1 )G B (ki,^)| 2 

qi 

+ ^2 E |M x (k,q 1 )G s (k 1 , / ,«)| 2 |M !/ (k 1 ,q 2 )G A (k 2) ^ 2 ])| 2 
qi,q2 

+ ^3 E |M x (k,q 1 )G B (ki,/^)| 2 |M„(k 1 ,q 2 )GA(k 2 ,^)| 2 

qi,q2,q3 

|M x (k 2 ,q 3 )G B (k3,ML 3 l)| 2 + - , (8) 

of the wavefunction on the A sublattice given by Eq. (4a) , and similarly N^b for 
the B sublattice. Naturally, N^l = 1 which yields the following suggestion: If 
the sum of the norms of the m first terms of the wavefunction fulfills the equation 
Y^n=o^kL — 1, then all terms with n > m in the Rciter's wavefunction can be 
neglected. 

In order to deduce what is the value of m for different values of the orbiton 

(n) 

energy J we calculated numerically N^ L for n = 0, 1, 2, 3 as a function of J on a 
mesh of 16 x 16 k points t, cf. Fig. 1. The obtained results do not depend on the 
sublattice index L. Fig. 1(a) shows the results for the orbital model Eq. (2a)- 
(2b) whereas Fig. 1(b) shows the results for the orbital model Eq. (2a-2c). We 
conclude that: (i) for both cases for J > OAt the norm of the wavefunction is close 
to 1 already for m = 3 and hence there are just up to three orbitons participating 
in the formation of a polaron, (ii) it seems that only few more orbitons are excited 
for smaller J, and (Hi) for the case with three-site terms depends slightly 

on the momentum k — it follows the dependence of the QP spectral weight 

a L (k) on k, with a maximum at k = (7r/2,7r/2), i.e. at the minimum 



of e(k), cf. Eq. (3), and (iv) the dependence of on k found only when the 
three-site terms are included demonstrates that the three-site terms are indeed 
responsible for the formation of mobile polarons. 
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t Obviously, in this procedure the QP spectral weight \a L | is calculated directly from the 
Green's function [6] and not from the normalization factor to the wavefunction. 
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Fig. 1. The norm N^j as a function of J for n = 0,1,2,3, from bottom to top re- 
spectively, as obtained for the orbital model (a) without and (b) with three-site terms 
Eq. (2c) included in the Hamiltonian. The results on panel (a) do not depend on k 
whereas dotted (dashed) lines on panel (b) show results for k = (0, 0) [k = (7r/2, vt/2)], 
respectively. All of the results do not depend on the sublattice index L. 

Summary — As a summary, it is instructive to compare the results of the 
realistic orbital t-J model (i.e. with three-site terms included) with those obtained 
for the spin Ising model and for the spin SU(2) symmetric model in Refs. [6]- [7]. 
On one hand, for all n the dependence of the norm on k resembles to some 
extent the spin SU(2) model. On the other hand, a closer look reveals that even 
for k = (0, 0) the dependence of on J is a concave function for all n just 
as for the spin Ising case and unlike in the SU(2) case where it can be a convex 
function of J for some n [7]. Hence, the detailed study of the Reiter's vavefunction 
of the orbital polaron confirms that the hole doped into the t 2g AO forms a mobile 
polaron just like the mobile spin polaron in the AF plane of high-T c cuprates [2] 
but its properties are truly distinct, in agreement with discussion in Ref. [1]. 
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